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In order to provide a formally correct thermodynamical description of inhomogeneous fluids valid 
on all length scales down to the classical limit we postulate that all extensive quantities have locally 
extensive analogues. We derive local thermodynamic relations for the canonical and isothermal- 
isobaric ensembles without limitations with respect to the nature of intermolecular interactions. 
The inhomogeneity of the system is shown to result in two distinct entropic contributions. The 
intrinsic contribution is due the specific distribution of local environments experienced by a single 
molecule while the extrinsic one originates exclusively from the local variation of the density from its 
bulk value. As an example the theory is applied to the surface tension of a spherical liquid droplet. 



The study of spatially inhomogeneous fluids is im- 
portant for a vast variety of applications, ranging from 
macroscopic surface and interfacial phenomena |l|, spa- 
tially confined systems such as fluids confined to channels 
of porous materials [2J , down to molecular solvation phe- 
nomena |3| . A proper theoretical framework for address- 
ing these problems should be capable to unambiguously 
describe local thermodynamical properties. The term 'lo- 
cal' is to be understood in terms of a function that is well 
defined at a certain point in space and in general also de- 
pends on the state of the system in the neighborhood of 
the point. The size of the neighborhood as well as the 
spatial scale, on which the description is valid, depend 
critically on the level of theory. 

Various theoretical attempts to describe 'local' ther- 
modynamic properties of inhomogeneous fluids have been 
proposed (traditionally called point thermodynamics) re- 
lying either on Kirkwoods integral equation [J] , spatially- 
dependent cluster expansions |5| or a density functional 
treatment [f| . In the fist two approaches [J, Q the defini- 
tion of local thermodynamic potentials necessitates pos- 
tulating the existence of a local pressure, which is jus- 
tified on spatial scales, on which local properties vary 
slowly, but it fails on smaller scales, on which one has to 
introduce a pressure tensor due to spatial inhomogeneity 
[6|. A proper description which enables a full coupling 
of local thermal and mechanical properties was made us- 
ing the density functional formalism [6] . While it consis- 
tently describes fluids near interfaces or confined in pores 
and also enables a direct extension to the description of 
transport of dense inhomogeneous fluids it builds on the 
assumption that the correlation function is short ranged 
Q • Thus it cannot be applied to describe systems on the 
scale of the correlation length or even below, such as of 
immediate interest, for example, in solvation phenomena 
[3| or molecules confined to nanotubes [8]. In these sys- 
tems one is often specifically interested in entropic effects 
(i.e. the constraining of molecular degrees of freedom). 
Inevitably one has to resort to some sort of correlation 
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expansion Q, which becomes impractical for inhomoge- 
neous systems and is therefore usually ad hoc truncated 
at the pair correlation level. As significant higher order 
correlations were found already in homogeneous simple 
liquids [10( 1 such a truncation can not be appropriate in 
the case of hydration [Tljj . More importantly, the use of 
thermodynamic concepts in the case when the dimension 
of the studied (sub)system is comparable to or smaller 
than the correlation length, is formally not justified. The 
immediate technological interest in such systems such as 
nanofluidic applications [12[ and wetting phenomena in 
nanotubes [13[, poses a pressing need for a general for- 
mulation of local thermodynamics on such small length 
scales. At present, such a formal framework is still elusive 
and even the question of mere validity of thermodynam- 
ical concepts on spatial scales comparable to or smaller 
than the correlation length is still under debate. 

Here we present an exact rigorous formulation of local 
thermodynamics of inhomogeneous fluids at equilibrium, 
which holds on arbitrary spatial scales down to the classi- 
cal limit without any limitations regarding intermolecu- 
lar interactions or the origin of the inhomogeneity. Based 
on the test particle picture we show that the specifica- 
tion of the order of correlations can be entirely avoided 
while they arc inherently built into the theory. Defin- 
ing locally extensive properties we derive a complete set 
of local thermodynamic relations in the canonical and 
isothermal-isobaric ensemble. An exact expression for 
the surface tension of a spherical liquid droplet is de- 
rived as a 'proof of concept' example. Finally, we discuss 
the physical meaning of local thermodynamic potentials 
on scales smaller than the correlation length. 

While the quantities pertaining to properties of single 
molecules (e.g. local density, (p(r))) have an unambigu- 
ous local definition, functions which are averages over 
the properties of larger groups of molecules (e.g. free 
energy density, f(r), or entropy density, s(r)), cannot 
be unambiguously defined but arc required to fulfill cer- 
tain constraints, typically the invariance of some of the 
lower moments of the functions [lj. We postulate that 
all nominally extensive quantities have well defined local 
analogues and therefore introduce the notion of a locally 



extensive quantity a(r) as: 

a(r) = {p(r))A\r), 



(1) 



where the superscript 1 denotes the property per single 
molecule located in the volume element dr centered in r. 
We describe the 'per molecule' or 'quasi-single' particle 
property within the test particle-insertion picture, where 
a local quantity is typically a function of the insertion 
energy of the virtual test molecule and assign the entire 
potential energy change to the test molecule. We show 
that such a choice leads to invariant values of integrals of 
local thermodynamic functions over space. We focus on 
excess thermodynamic properties since the corresponding 
ideal contributions are trivial to evaluate. 

In the canonical ensemble we begin by defining f e x(r) 
in terms of the excess free energy per molecule located 
in the differential volume element, F ex (r), which can be 
understood as the change in excess free energy upon in- 
sertion of an additional molecule at position r at con- 
stant T and V integrated over all orientations of the in- 
serted molecule and all configurations (positions as well 

as orientations) of the remaining N molecules. F ex (r) 
is expressed in terms of the corresponding configura- 
tional integrals of the (N + 1)- and iV-particle systems: 

T ex (r) = -fc B rin ZjY +^ (r) . We denote the intermolecu- 
lar interaction energy of the N molecules as Un and the 
interaction energy of these molecules with the external 
field (the potential due to the presence of the solvent) 
as ipN- The potential energy of the inserted (N + l)-st 
molecule, including both the interactions with the re- 
maining N molecules and the interaction with the exter- 
nal field, is denoted as ip. Thus, the potential energies 
read Un + 4>n and ip + U N + 4>n = U N +i + "4>n+i in 
the N and (N + 1) system, respectively. The configu- 
ration of the z-th water molecule is specified by stating 
its position and orientation, s^ = (rj,f2j). Written out 
explicitly, f e x(f) reads: 



f ex (r) = -p- 1 (p(r))ln( / dfl exp(-^(r, fl)) 



N 

(2) 

where /3 = 1/fcgT and the subscript N denotes the aver- 
age over all configurations of the remaining N molecules. 
In the case where the system is homogeneous (de- 
noted with index 0) Eq. (jSJl reduces to the bulk den- 
sity multiplied by F ex0 , the latter term being equal 
to Widom's formula for the excess chemical potential 

PbF £x0 = psPex- In the case where the external po- 
tential, ip(r), is spherically symmetric, Eq. © splits 
into an external field contribution, (p(r))ip(r), and an 
intrinsic term, — /3^ 1 (p(r)) In ^JdOexp(— f3u(r, fi))) N , 
with ip — ipN+i — VVj and u = Un+i — Un- This 
separation establishes a connection with density func- 
tional theory (DFT) and the Yvon-Born-Green hierar- 
chy. Namely, it is known that the chemical potential of a 
inhomogeneous system splits into the intrinsic chemical 
potential, Pi n t, and the external field contribution, i.e. 



P = Pint{r) + ip{r). Explicitly, pi n t{r) is the functional 
derivative of the intrinsic free energy functional [l|, Il4| , 

s(p(r\\ — (J-intir). The intrinsic free energy functional 
can be divided into ideal and excess parts and the latter 
can be identified formally as the excess intrinsic chemi- 
cal potential, M?nt( r ) — g/"}^)) ( see f° r example fllll4l|) 
or equivalently the singlet direct correlation function, 
c x (r) in integral equation theories and the Yvon-Born- 
Green hierarchy [l|, [lj, [l5| . Note that if instead we had 
ip = ip(r, fl) a DFT free energy functional would read 

F ex = J dr J dfl(p(r, tt)) N F ex (r, $1). A factorization of 
the chemical potential into intrinsic and external contri- 
butions as to be used in a DFT-like approach would then 
be rather useless. Meanwhile a generalization of the first 
member of the YBG hierarchy can in fact be obtained 
formally by using the invariance of the system under si- 
multaneous rotation of the molecules and the coordinate 
system [lfjj. 

The excess entropy density, s ex , is obtained from the 
thermodynamic identity S ex — — gj^- Taking the par- 
tial derivative inside the integral and using the subscript 
f2 and a shorthand notation for the angular integration, 
(•••)n = J ■■•dfl we obtain Sex = - J -§ffex{r)dr = 
J s ex {r)dr or explicitly 

sex(r) = i (/T^r)) ln(exp(-)Mr, n ))n)jv 

d 



- (p(r))—ki(eM-Mr,n))n)N)- (3) 

The first and third terms in Eq. ^ represent the con- 
tribution of intrinsic excess entropy of single molecules, 

S ex (r) = Jt i and expresses the number of de- 
grees of freedom of single molecules at r, while the 
second term can be understood as an extrinsic excess 
entropy density, —F ex (r) ^ , and is solely due to 
the fact that the system is inhomogeneous (it is zero 
for a homogeneous system). Due to the indistinguisha- 
bility of molecules the local density can be written as 
(p{r)) = (2j=i ${ r i — r )) — N{S(ri —r)), and we obtain 
after straightforward algebra that 



^(p(r)) - (p(r))(W N ) N - N(WNS( ri - r)) N - (4) 

with Wn = Un + ipN ■ Eq. ^ gives a vanishing change 
in particle number when integrated over the volume as 
it should (since we work in the (N, V, T) ensemble) . By 
defining the average in the TV-particle ensemble where 
molecule 1 is fixed at r (or equivalently, an average over 
all configurations in which one of the molecules is located 
at r) as 



(O) 



(OS(n-r)) N _ (Oexp(-P<p{r,tt))n)N- 



N 



(5(ri-r)) N ~~~ (exp(-Pip(r,fl))a) N - 



(5) 



we find that the following relation holds 



J9_ 
dp 



(p(r)) = (p(r))((W N ) N -(W N r N ), 



(6) 



Eq. (j6|) expresses the fact that the change of local density 
with temperature depends on the difference between the 
potential energy of the system and the average energy of 
configurations where one molecule is fixed at r. Eq. (JOJ) 
tells us that the local density will decrease with increas- 
ing temperature if the average energy of configurations 
with one molecule fixed at r is more favorable than the 
potential energy of the system. It indicates that a lo- 
cally increased density with respect to the bulk density 
(i.e. the density sufficiently far away from the source of 
the perturbation) is disfavored in terms of the excess en- 
tropy as it naturally decreases the configurational volume 
in phase space. Elementary algebraic manipulations give 
the following expression for the intrinsic excess entropy 
per molecule 



-rjintr,! 



(r) = fc B ln(exp(-j8<p(r,n))n)jv+- 



(W)n+i - (W) 



N 



(7) 
where we have dropped the index for the potential energy 
indicating the number of molecules involved as it is equal 
to the one denoting the ensemble average. Using Eqns.© 
and (J7J Eq. ([3]) can be rewritten in the final form 



Sex(r) =7p(p(r)) { ln(exp(-/fy(r, il)) n ) N x 

{fr 1 + (W) r N -(W) N } + 
(Wy N+1 - (W) N ) 



(8) 

Here, using E = F + TS, the expression for the excess 
energy density, e ex (r), can be trivially written as: 

e ex (r) = <p(r)>(ln(fflq>(-Mr,n))oM(W% ~ < W M 

+ (W) r N+1 -{W) N ). (9) 



Eqns. ((2|), dHJ and ((HJ represent the fundamental lo- 
cal equations for the canonical ensemble. We see that 
the nominally extensive property exactly generalizes to 
the locally extensive property for inhomogeneous sys- 
tems. The local equations contain an intrinsic ('per single 
molecule') part and an additional extrinsic contribution, 
(p(r)) \n{exp(-/3ru(r, Q.)) n ) N ( w )n~(w)n ^ due tQ the h> 
homogeneity of the system. The extrinsic component 
expresses the energy redistribution. It is caused by a 
local density change with respect to a homogeneous sys- 
tem due to the redistribution of molecular positions at 
given local density of energy states per molecule. The 
intrinsic part reflects the local density of energy states 
per molecule and the change of potential energy of the 
system when a molecule is added to the ensemble at the 
specified position. We deliberately avoid defining the lo- 
cal pV term, as the scale-independent definition of the 
local work is ad hoc not possible because of its inherent 



non-local nature. However, we show in the following that 
the appropriate local pV analogue arises naturally in the 
(N,p,T) ensemble. 

Analogous to the (N, V, T) ensemble we again assume 
that the Gibbs free energy density is locally extensive, 
g ex (r) = G ex (r)(p(r)), where G ex (r) is the excess Gibbs 
free energy per molecule located in the differential vol- 
ume element and (p(r)) — (J2i=i <K r — r i)) lii tne local 
density (at constant N, p and T) in the presence of an ex- 
ternal field tp. The total Gibbs free energy per molecule 
located in the differential volume element centered at 
r is defined in terms of the configurational isothermal- 
isobaric partition functions of the N and N + 1 systems: 

Gl(r) = -fc 5 rin Ajv A ^ (r) [Hi. The excess Gibbs free 
energy density can be written as, 



g ex (r) = -r 1 (p(r))ln 



(V J dJ7exp(-/3^(r,rj))) 



/v 



(V) 



N 



(10) 

In the limit of small volume fluctuations, AV/V — > 0, the 
volume terms cancel each other and Eq. fTU)) reduces to 
Eq.([2]). The fluctuating volume term will be important 
where volume fluctuations are expected to be large. Im- 
portant examples include dewetting phenomena of meso- 
and macro-scopic surfaces [18], [l9| , hydrophobicity on the 
mesoscale [2(| and dewetting transitions in nanochannels 

m- 

The factorization of g e x(r) into excess enthalpy den- 
sity, h ex (v), and s ex (r) can be achieved by introducing 
the average volume (V)^f + i in a (N + l)-particle ensem- 
ble in which the (N + l)-st molecule position is fixed at 
r using Eq. (|3J|: 



9ex(r) 



/3- 1 (p(r))ln(cxp(-^(r))o)Ar 
(V)n+i 



/?->(*■)) In. 



(V) 



N 



(11) 



where the first term is recognized as the (Hclmholtz) free 
energy density and the second term reflects the relative 
response of the system volume to the insertion of an ad- 
ditional molecule at a point r given fixed T and p. Com- 
paring Eq. (ITT]) with the relation between the Gibbs 
and Hclmholtz free energies we recognize that the sec- 
ond term in Eq. (fTTj) is the exact local analogue of the 
pV term and immediately demonstrates the fact that the 
work is intrinsically performed non-locally. At the same 
time it also demonstrates that the total work has a well 
defined local contribution, which is exactly locally exten- 
sive. h ex (r) is obtained directly as g e x(f) +Ts ex (r) and 
all three together constitute the exact fundamental local 
relation for the isothermal-isobaric ensemble. 

Meanwhile, the chemical potential is strictly constant 
throughout the system and can only artificially be given 



a local physical interpretation, p e 



I P l e°x c (r)dr, 



terms of an average contribution per molecule p ex — 
N- 1 J g ex (r)dr. 

The locally-extensive formulation of statistical thermo- 
dynamics reduces to the standard description in case of 



a homogeneous system. On the other hand it must also 
reduce to the classical bulk thermodynamic formulation 
upon spatial integration. By noting that the Helmholtz 
free energy is related to p via F = Np—pV = N<fi, where 
<j> is the (Helmholtz) free energy per molecule. Recalling 
that the ideal gas contribution to F in the inhomoge- 
neous local picture is F^ = (3~ 1 J(p(r))(\n(p(r)) — l)dr 
we obtain for F = F^ + F ex , disregarding kinetic terms, 
/3 _1 f{p(r)) ln[(p(r))/(exp(— (3ip(r, n)))jvjdr. According 
to the potential distribution theorem [22| the logarith- 
mic term is constant, leading to F = N<j>, as expected. 
However, this suggests that Widom's definition of the 
chemical potential for an inhomogeneous system (within 
the framework of the potential distribution theory) [22| is 
in fact the free energy per particle. Clearly, for N — > oo, 
p and 4> become equivalent. 

As a case study we present an exact derivation of the 
surface tension of a single component liquid droplet in 
absence of external fields such that the liquid/gas intre- 
face is on average spherical. Addressing this problem in 
the framework of statistical mechanics using either me- 
chanical or point thermodynamical approaches is known 
to introduce a degree of arbitrariness into the problem 
such that the surface of tension (the surface in which the 
tension acts) is said to be ill-defined [l(. For notational 
convenience we assume that intermolecular interactions 
are spherically symmetrical but may include many-body 
contributions of arbitrary order. The system contains 
internal liquid (/) and external gaseous (g) bulk phases 
separated by a non- uniform interfacial region (s). Choos- 
ing the radius of the dividing surface, r s , determines the 
volumes of the bulk phases V = Vi + V g . Constancy of 
total number of molecules demands N = Ni + N g + N s — 
An J r 2 (p(r))dr. Only excess quantities contribute to the 
surface tension and F ex of the system can be written in 
terms of bulk and surface contributions F = Fi+F g + F s . 
Depending on r s the values of N s and F s can formally be 
positive or negative. By definition the surface tension, 
a, represents the excess free energy of the dividing sur- 
face per unit of its area: a = (A-Kr 2 s )~ l (F — (F g + F{)). 



Extensivity, Fi 



Kg) 



N l(g) F ex,U 



(where the subscript 
denotes the homogeneous iV-molecule ensemble of bulk 
phase 1(g) at the same T and V), allows us to write 



a = r„ 



dr- 



r 2 ( P {r)){F 1 ex {r)~F 1 f x ^ 
r 2 (p(r))(Fl(r)-Ft 9 fi )dr}. (12) 



We have a — o~(r s ) and it is known that at the actual 
surface of tension the formal derivative vanishes at rt, 



da/dr s \ rs=rt = [1]. Applying the chain rule and making 
use of the second fundamental theorem of calculus we can 
perform the derivation and obtain an implicit equation 
for 7V 



-h 1 ^ -fiUK 



■(p(n))(F e :~F e ' x ) = 



r 2 {p(r))$ l (r)dr 



r 2 (p(r))§ 9 (r)Ar, (13) 



Once we know (p(r)), 



where fl'W (r) = In ^'^7^ 

Fex,o an d ( ex P( — Pf( r )))N, all of which can be read- 
ily calculated from Monte Carlo or Molecular Dynamics 
simulations, we can compute the surface tension as 



-(p(r)) {d l (r)H{r t -r) + &(r)H(r - n)} dr, 

r t 

where H(x) denotes the Heaviside function. Analogous 
equations are obtained in the case of a flat interface, 
where the surface of tension is chosen such as to satisfy 
N s = 0. Eqs. ([T2"T) and (fT4| show that the local density of 
surface free energy, (local density of a A), is also exactly 
locally extensive. The above derivation can be viewed as 
the exact accomplishment of the van der Waals original 
approximate program for flat interfaces based on point- 
thermodynamics [lj, |23| . The present method is free of 
arbitrariness and avoids the use of local pressure tensors, 
which are known to be non- uniquely defined |lj. 

In this letter we propose an exact general formula- 
tion of local thermodynamics of inhomogeneous fluids in 
the canonical and isothermal-isobaric ensembles, which is 
valid on all length scales without limitations with respect 
to nature intermolecular interactions. By introducing the 
notion of a locally extensive quantity we have shown that 
the thermodynamic potentials, which are inherently non- 
local, all have a well defined local physical meaning. Our 
results show that there are two distinct physical conse- 
quences of inhomogeneity, an intrinsic one due to the 
specific properties of the local environment giving rise to 
a specific local density of states of a molecule, and an ex- 
trinsic one, which is explicitly due to the inhomogeneous 
distribution of molecular positions throughout the sys- 
tem. The theory was used to derive an exact expression 
for the surface tension of a spherical liquid droplet. Our 
results pose the formal basis for a general local thermody- 
namic description of inhomogeneous fluids ranging from 
macroscopic interfacial phenomena, molecular solvation 
to such as encountered in modern nanofluidic devices. 
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